CORRELATION BETWEEN POLE LOCATION AND 
ASYMPTOTIC BEHAVIOR FOR PAINLEVE I SOLUTIONS 



OVIDIU COSTIN 

Abstract. We extend the technique of asymptotic series matching to exponen- 
tial asymptotics expansions (transseries) and show that the extension provides a 
method of finding singularities of solutions of nonlinear differential equations, us- 
ing asymptotic information. This transasymptotic matching method is applied to 
Painlcvc's first equation. The solutions of PI that are bounded in some direction to- 
wards infinity can be expressed as series of functions, obtained by generalized Borcl 
summation of formal transseries solutions; the series converge in a neighborhood of 
infinity. We prove (under certain restrictions) that the boundary of the region of 
convergence contains actual poles of the associated solution. As a consequence, the 
position of these exterior poles is derived from asymptotic data. In particular, we 
prove that the location of the outermost pole x p (C) on R + of a solution is mono- 
tonic in a parameter C describing its asymptotics on antistokes lines 1 , and obtain 
rigorous bounds for x p (C). We also derive the behavior of x p {C) for large CeC. 
The appendix gives a detailed classical proof that the only singularities of solutions 
of PI are poles. 



1. Introduction and notations 
We consider the solutions of Painleve's first equation, PI 

(1.1) y" = Qy 2 + z 

and study the connection between their asymptotic behavior and the location of their 
poles. Classical asymptotic analysis and isomonodromy theory [1] have not yielded 
so far this type of information. 

The analysis and proofs in the present paper are mostly self-contained, but the 
method that we use, as well as the nature of results obtained, are much easier to 
interpret from the more general perspective of exponential asymptotics. We start by 
briefly presenting some concepts and results in this area which are relevant to the 
problem under discussion. 

Jean Ecalle has introduced and developed the theory of analyzable functions and 
tools of generalized resummation of transseries (generalized exponential-power series) 
as a natural and very comprehensive generalization of the theory of analytic func- 
tions [2, 3, 4]. The grand picture is the construction of a class of functions, closed 
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under all imaginable algebraic and analytic operations, which are into an isomorphic 
correspondence with the space of transseries. In a sense, the space of transseries 
represents the formal closure of power series under all operations. A wide class of 
analytic problems admit formal solutions in the space of transseries (obtained in a 
more or less straightforward way). Transseries solutions contain ample analytic infor- 
mation and in many cases (for instance generic linear or nonlinear differential systems 
[8]) the actual solutions of the problem can be reconstructed from their transseries 
by modern resummation methods. 

1.1. Notes on asymptotics beyond all orders for PI. Proofs of the various 
technical statements in this section can be found in [8], in the more general context 
of analytic rank-one differential systems. 
The change of variables 

(1.2) iVQz^y{z) = -1 + -i- - h(x); 30x = (-2Azf /A 

25x 2 

brings (1.1) to the (normal) form 

(1.3) h " + L h '- h -^-^L = 
y ' x 2 625x 4 

There is a one-parameter family of solutions of (1.3) that tend to zero as x — ■> +oo, 
i6K, and they have a common asymptotic series (cf. also § 3): 



(1.4) h(x) ~ ho(x) = ^c k x k 

k=A 

The classical expansion (1.4) does not distinguish different solutions in this family. 
However there exists a one complex parameter family of transseries solutions, of which 
h is only a special instance: 



(1.5) h(x, C) = ho(x) + C k x- k ' 2 e- kx ~h k {x) (C e C) 

k=i 

In (1.5), hk are formal power series, hk = Ylm=o c kmX~ m , that can be computed 
algorithmically by substituting (1.5) into (1.3) and identifying the coefficients. With 
the choice c w = 1, which is possible for (1.3), the series hk do not depend on C. Every 
hk is factorially divergent. There is nevertheless a very close connection between (1.5) 
and the family of actual solutions of (1.3) which decrease on R for x — > +oo: any 
such solution can be reconstructed by generalized Borel resummation of h(x, C) for 
some C. Denote by CB the generalized Borel resummation operator (as defined in 
[8]). There exists xo > (independent of k) so that the functions hk = CBhk are 
analytic in the half plane H XQ = {x : 9?(x) > x }. In addition, hk ~ hk as x — > e l9 oo 
in H Xo , 9 E (— 7r/2, 7r/2). For any C 6 C, the function series 
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(1.6) h(x; C) = £Bh (x) + C k x~ k,2 e- kx CBh k {x) 

k=l 

converges for $l(x) > X\ = X\(C) to a solution of (1.3); conversely, any solution h 
which vanishes for real x — > +oo is necessarily of the form (1.6) for a unique C. 
LB commutes with the algebraic operations, complex conjugation (thus a series with 
real coefficients sums to a real function), with differentiation, and has good continuity 
properties. Along directions in C other than IR + similar results hold. 

In addition, the constant C in (1.6) is related by a simple formula to the large 
\x\ behavior of h(x;C) on the antistokes lines ±iR + (cf. (2.8)). We show that for 
any C > 0, the boundary of the domain where the series (1.6) converges contains an 
actual pole of h(x; C), much as the boundary of the convergence disk of the Taylor 
series of an analytic function must contain a singularity. 

1.2. Heuristic discussion of transasymptotic matching. The value of C is re- 
lated to the position of the outermost singularities of solutions of PI. The following 
heuristic argument illustrates the reason. 

For simplicity, we assume that C is large._ When x is very large, the exponentials in 
(1.5) are smaller than the leading term in ho, and the behavior of h(x; C) is roughly 
c 04 x' 4 . 

For smaller values of x, such that 1 ^> \Ce~ x x~ x l 2 \ ^> the relative magnitude 

of the terms of (1.5) changes: the largest term becomes Ce~ x x~ 1 ^ 2 (instead of co^x" 4 ), 
the next order of magnitude is the term C 2 e~ 2x x~ 1 C2o (not Ce~ x x~ 3 / 2 Cu) and so on. 
In this intermediate regime, the terms of (1.5) need to be reshuffled accordingly, and 
h(x; C) is represented by 

(1.7) h{x; c) = jr — (jr , Xxl k /y C)k ) ■= £ m( x m 1 ) 

m=0 \k=0 v ' ' / m=0 

It turns out that, for x not too small, the inner sums converge (while, obviously, 
the expansion ft as a whole is still divergent), and G m {t) are analytic. In fact, 
for (1.3), G m (t) are rational functions and Go has a particularly simple expression, 
Go(t) = t(t — 1/12)~ 2 (see the proof of Lemma 8). Moreover, for any m > the 
function G m (t) has only one pole in the right half plane, at t — 1/12. 

The representation (1.7) is normally accurate as long as no Stokes phenomena 
occur and the following asymptoticity condition holds 

(1.8) \G m+1 (e x x 1 / 2 /C)x- m - 1 \ < \G m (e x x l / 2 /C)x~ m \ 

A formal calculation shows that (1.8) holds for x in a region 1Z that starts well to 
the left of the array of poles of g(x) = G (e x x 1 ^ 2 /C), located at xn, where 



(1.9) 



x N + - \nx N = ln(C/12) + 2Niri; N E Z 
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and extends all the way (small neighborhoods of the poles themselves excluded) to 
di(x) — > oo. The region 1Z is thus substantially wider than the region where (1.4) is 
valid. 

Accepting for the moment that (1.7) is the asymptotic expansion of some solution 
h in 71, (1.7) can be used to estimate §(j h) on a circuit around a pole xn- The 
calculation gives <f(f h) = —12ni + 0(l/ InC) for large |C|, showing that h is singular 
at some point x p (C] N) = x^+o(l). x p (C] N) can be estimated to all orders in 1/lnC 
from (1.7). 

On the other hand, if x is large, and away from xn, h(x; C) must be bounded (by 
(1.7)). The location of the poles of h in 1Z should thus be x p (C; N) = xn + o(l). 

This method gives information not only about the location, but also about the 
nature of the singularities of h(x; C) in 1Z. A straightforward calculation based on 
the Cauchy integral argument explained above indicates, to any order in 1/ InC, that 
x p (C; N) are double poles. 

The fact that PI is an integrable equation plays no role, except that G m can be 
computed in closed form. 

We note that in the matching process all the terms of the transseries (1.5) are nec- 
essary, so that the matching is transasymptotic. The procedure makes the connection 
between the behavior of h in the pole-free region near infinity, and its singularities 
for large x in C. 

In the next sections these conclusions are made rigorous, and further results are 
proved on the connection between C and the position of x p (C; 0), for C not necessarily 
large. The proofs yield global estimates for (1.6) and (1.7). 

We can also interpret (1.7) from the perspective of multiscale analysis. Apparently 
an expansion of the form c m (x)x~~ m (where the coefficient of x~ m may depend 
on x) is subject to substantial arbitrariness. Not so, however, if c m (x) are required 
to be sufficiently regular and the expansion to be valid over a wide enough range 
of x. This is a crucial point in multiscale analysis. We note that (1.5) determines 
the asymptotic scales x~ l and Ce~ x x l l 2 . Once the scales are found, the terms in 
(1.7) are more conveniently determined by inserting (1.7) in (1.3) and identifying the 
coefficients of x~ m . 

The proper choice of the asymptotic scales is essential. For instance looking for a 
representation of the more general form h a = J2x~ m G m (Ce~ x x a ), only for a = 1/2 
will the expansion be asymptotic (even formally) for large x. For any other value of 
a, G m (Ce~ x x a ) = 0(x m ) instead of 0(1), and then h a is rather meaningless. 

2. Main results 

Notations. By Lemma 2 in § 3 there is a one (real) parameter family of solutions of 
(1.3) which are real- valued for x > const, on M + and such that ho(x) — > as x — > +oo. 
Their behavior for large x is ho ~ c^x~ A + C6X _6 (1 + o(l)) where C4 = —392/625 and 
cq < 0. We choose ho to be one such solution. (The choice h = CBho has however 
important additional properties). We let A > 9/4 be such that, for x > A, 



(2.1) 



^o(^) ^ ( — x 4 , c 4 x 4 ) 
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Proposition 1. (i) Any solution h of (1.3) such that \h(x)\ < 1 for large positive x 
is given by 



(2.2) h(x; C) = C k x- k ' 2 e- kx h k (x) 

k=0 

for some CgC, where the functions h k (x) are determined recursively as solutions of 
linear differential equations (cf. (3.2)) and do not depend on C. 

The series (2.2) converges for all x > xq{C) = m&x{A,xi} where x\ + |lnxi = 
ln(|C|/12). 

(ii) The functions hk satisfy the inequalities 

k 

(2.3) < h k (x) < lim h k (x) = (k e N and x > A) 

x—>oo 12 

Furthermore, 

lWh k (x)' 1/k = 12 + 0(x~ v ) (as x -> +oo) 

(Hi) (Outermost pole of h(x; C) on R+.) If C > C = 12min t > A y/te^l - J)" 1 , then 
h(x; C) has a double pole at the point x p = x p (C; 0), the unique solution of 

(2.4) i'/VIEA^i)- 1 /' = C 

for x > A. The function h(x; C) has no other singularities on [x p , oo). 

(iv) (Monotonicity and bounds for x p .) For C > C , x p (C; 0) is increasing in C and 
satisfies the inequalities 

1 1 ( 9 \ 

(2.5) x p + -lnxp < ln(C/12) < x p + -\nx p — In I 1 — - — 

2 2 \ 4:X p J 

Thus, x p has the asymptotic behavior, 

(2.6) x p = ln(C/12) - iln(ln(C/12)) + 0(1/ In C) (C -> +oo) 

(v) (Poles of h(x; C) in C, for large C complex. ) Let N e N and e > 0. As \C\ — > oo ; 
h(x; C) has double poles at the points: 

(2.7) x p (C; k) — ln(C/12) - ^ ln(ln(C/12)) + 2km + o(l) (k eZC\ [-N, N]) 

and no poles in the region {x : \$s(x)\ < N, > 9?(x p (C; 0)) + e}. 

(vi) (Connection to behavior on antistokes lines. ) For the special choice ho = CBh 
we have h k = CBh k , and 



(2.8) 



h(x; C) ~ (C ± ^aT^e"* osi^ ±ioo 
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where C G C is the same as in (2.2) and S is the Stokes constant of the equation 
(1.3). 

Notes. (1) The property that the solutions which vanish in some direction can be 
expressed, as in (2.2), as a power series in the constant C of (2.8) is very general, cf. 
[8] , and explains why C is also a natural parameter for the location of exterior poles. 
Integrability of PI plays no role here. A formula of type (2.4) is expected to be quite 
general and is not a connection formula stricto sensu, because it involves limits. For 
PI however, there are indications that these limits can be calculated in closed form. 

(2) For the choice h = CBh it can be shown, using the methods in [9], that (2.1) 
holds for x > 4/3. 



3. Proofs and further results 

Before proving the main proposition we describe some asymptotic properties of the 
solutions of (1.3) which are bounded as x — > +oo. 

Lemma 2. (i) If h is a solution of (1.3) and \h{x)\ < 1 for large positive x, then 

oo 

(3.1) h~h = J2 c >< x ~ k ( x eie °° with 9 e 

fc=4 

(the first few terms of h are -§§a;- 4 - ^x" 6 - ""oS^ " 8 -■■■)• 
(ii) There is a one (complex/real) parameter of such solutions (complex/real valued 
on if hi and hi satisfy (3.1) then hi — hi — Ce~ x x~ 1 ^ 2 + o(e~ x x~ l l 2 ), for some 
C G C, as x — > oo. If in addition C = then hi — hi. 

Proof of Lemma 2. (i) is standard (see, e.g. [6]); it is also easy to reprove in our 
case, since (1.3) can be rewritten as 

h" + -h'-(l + ^]h= -h 2 ~-^h + -^- A := g(x) 
x \ Ax 2 ) 2 4a; 2 625x 4 yv ; 

which, in integral form 

h(x) =x- 1 / 2 (e x f e-H^ 2 g(t)dt + e- x f e^t 1 ' 2 g{t)ti\ 

\ J a J oo J 

is a contraction in the unit ball in the sup norm, for large x. If a G R then h is 
real-valued. 

Part (ii) is proven without difficulty in essentially the same way, by writing an 
integral contractive equation for 5 — h 2 — hi (see [8] for details under more general 
assumptions). □ 
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3.1. Properties of (2.2) on M + . Let h and A be as in § 2. Formal substitution 
of (2.2) in (1.3) yields 



(3 



.2) El - (l + h (x) - JLW = y^e~ kx J2 h i h *-i : = R *&) ( k ^ !) 

0<j<k 



where H k (x) = x (fc 1)/2 e kx h k (x). 

Step 1. Study of the homogeneous part of (3.2) (same as (3.2) for k = 1): 

(3.3) y" - (l + h (x) - JL) y = 

Lemma 3. (i) Eq. (3.3) has two solutions y+,y~ such that: y±(x) > for x > A, 
U±( x ) — e ±x (l + 0(x~ 1 )) as x — > +oo and their Wronskian (y~y' + — y'_y+) is W = 2. 

If y(x) solves (3.3) and y(x) = 0{e~ x ) as x — > oo then y = const. y-. 
(ii) The function u(t, s) = y + (t)y-(s) — y-(t)y + (s) is positive for t > s > A. 

Proof, (i) Except for the explicit region of positivity of y±, this is standard asymp- 
totic theory (see e.g. [6]; note also that W = const and W — > 2 for x — > oo). We 
prove the statement only for the argument for y + is similar. Associating to (3.3) 
the integral equation 

(3.4) f(x) = 1 - jH e*~ s sinh( S - x) (± - h ^j f(s)ds 
where f(x) = e x y(x), we define 

(3.5) J(V)(x) = J e x - s smh(s-x)(^-h ^JV(s)ds 
on L°°(a;o, oo), for x > A. 

Remark 4. || J\\ L oo^ L oe < ( 1 H ^ J and J positivity preserving. Thus (3.4) 

8xq \ Sxq J 

has a unique solution f\ in L°°(A, oo); moreover, < /i < 1 and /i(x) ~ 1 + 

Er=i(- i ) j 'i^i^(^^+ o °)- 

Proof. The first part is straightforward. Now, ||/i||oo < (1 — ll^1l) -1 ||l||oo thus 
|/i| < | implying \f\ — 1| < |J(1) < |, in particular f\ > 0. Since J preserves 
positivity, we have f\ = \ — Jf\ < 1. The asymptotic expansion of /i follows from 
the representation f\ = (J2 < JL (—jy) 1- □ 

For the proof of Lemma 3 (ii) note that u(t, s) is a solution of n tt = (l + h (t) — \t~ 2 )u 
with initial condition u{s, s) = 0; u t (s, s) = 2. Since 1 + h (t) — \t~ 2 > for t > A, 
we have n(t, s) > for all t > s. □ 
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Lemma 5. There is a unique solution of (3.2) with Hi = y_ such that for k > 0, 
H k {x) = 0(e~ kx ) as x — > oo. Namely, the H k {x) are inductively given by 

i r°° 

(3.6) H k {x) = -j u(t,x)R k (t)dt, k>2 

In addition H k (x) > for x > A. 

Proof. Immediate induction using Lemma 3. Uniqueness follows from the behavior 
for large x of the solutions of the homogeneous equation, y± (Lemma 3). □ 

Step 2. Upper bounds for H k . We have, by (3.2) and Lemma 5 



(3.7) H k (x) = I sinh(t - x)R k (t)dt - / sinh(t - x) - h (t) ) H k {t)dt 



< 



/ sinh(t - x)e" kt R k (t)dt 

J X 



From Lemma 3 and Remark 4 we get Hi(x)e x G (0,1) on (A, oo). Assuming by 
induction that H m (x) < 12 ™-i ^~ (fc ' 1)/2 e~" 1:r if 1 < m < k, we get 



(3.8) H k {x)<\ r^{t-x)r k We- kt %T^< 

0<j<k 

(k 3 - k) (fc -i)/ 2 r ginh(t _ x)e -kt dt = k x -{k-i)/2 e - kx 
12 k J x v ; 12 fc - x 



□ 

Step 3. Monotonicity 

Lemma 6. The functions H k (x) and F k (x) = x~ 1 ^ 2 e~ kx h k (x) , k>l, are decreasing 
for x > A. 

Proof. The functions F k satisfy the equations 
(3.9) Ft = (1 + h )F k - ±F' k + l - FjFk-j (k > 1) 

0<j<k 

By (3.8) and Lemma 5, the F k are positive and vanish as x — ■> oo. On the other 
hand, F k can have no extremum for x > A (otherwise, by (3.9) all extrema would be 
minima), thus are decreasing in x. For H k the proof is very similar. □ 

Step 4- Lower bounds for H k . 
Lemma 7. For x > A we have 

(1 \ / Q \ ( m ~l)/ 2 
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From Remark 4 we get / = 1 - Jf > 1 - Jl > 1 - (1 - (8a;)- 1 )- 1 for x > 5/3. Thus 
H\(x) > (1 — (8a;)- 1 ). Since the functions H k are decreasing, we have 



dt sinh(t - x)R k (t) - H k (x) J dt sinh(t - x) I — - h 



(t) 



Thus 



(3.12) H k (x) > (l - jT sinh(t - x)R k (t)dt 

From this point on, the proof of Lemma 7 follows the same track as Step 2, noting 
that (1 — £ _1 )( fc_1 )/ 2 is increasing in t and that f°°sinh(t — x)t~ < - k ~ 1 ^ 2 e~ kt > (1 — 
x^ 1 )(k 2 - l)~ 1 x~ { - k ~ V) / 2 e~ kx . ' □ 

Step 5. Completing the proof of Proposition 1 (i) and (ii). With C 6 C given, 
(2.2) converges for x > xo(C), by (3.8). Furthermore, since the H k are decreasing, if 

(2.2) converges at x — x±, then it converges uniformly on [xi,oo) (together with its 
derivatives, in view of the equality in (3.7)) and thus is a solution of (1.3). By (3.8) 
and (3.10), we have h(x;C) — ho(x) = Ce~ x x~ 1 ^ 2 (l + o(l)) for large x. Therefore, 
by Lemma 2, if g is solution of (1.3) which vanishes as x — > oo, then g = h(x; C) for 
some C e C. □ 
Step 6. Divergence of the series (2.2) implies singularities of the function h(x; C). Let 
C > Co andxdiv '■= inf {x| Y^k=i C h F k (x) < oo}. We have just shown that (2.2) solves 

(1.3) for x > Xdiv By the definition of Co and (3.7) we have Xdiv > A. We now want 
to show that the solution h(x; C) is nonanalytic at xa v . Let Sn{x) = ^2 k=1 C k F k {x). 

Lemma 8. Assume h(x;C) is a C 2 (x 1 ,oo) solution of (1.3), and xi e (A,Xdi V ]- 

(i) The function Sjy(x) = h(x;C) — h — Sjy(x) is positive for all x > X\ and 

N e N+. 

(ii) xi = x div = x p (C; 0). 

Proof, (i) We have, by Proposition 1 (ii) 5^(x) = J2j >N C^ Fj(x) > for x > Xdiv 
Now, since the function 5n satisfies the equation 

(3.13) S'^ = \s 2 N + (l + S N + h )6 N -^- + \ J2 F ^ 

i,je(i,N-i) 

i+j>N 

we see, as in Lemma 6, that 5n is decreasing, therefore positive, for all x > x±. 

(ii) By (i), h(x;C) is unbounded at x^ v . But the only singularities of h(x;C) for 

x 7^ are double poles (see the appendix). □ 

Step 7. Concluding the proof of Proposition 1, (Hi) and (iv). The monotonicity of 
x p (C; 0) in C is now immediate, while (2.5) follows from Lemma 8, (3.8) and (3.10). 
□ 

Note. Formula (2.4) in itself is not a practical numerical way to calculate x p {C), 
given a particular C. An efficient (and rigorously controllable) method is to use the 
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procedure in this section, to obtain successively better upper and lower bounds for 
Hk- 

3.2. Correspondence between (1.6), (2.2) and (2.8). 

Lemma 9. If h = CBho then hk = CBhk for all k G N. Thus the representations 
(1.6) and (2.2) coincide. Formula (2.8) holds. 

Proof of Lemma 9 and of Proposition 1, (vi). We let gi i2 (x,C) = h(x;C) as given 
by (1.6) and (2.2), respectively. In [8], formula (2.8) for gi(x,C) and the estimate 
gi(x,C) — ho(x) = Ce~ x x~ x l 2 (\ + o(l)) are proved in a general setting. Lemma 2 
then implies, as in Step 5 § 3.1, that gi(x; C) = g2(x; C). □ 

Note. Using Lemma 9 it follows from [8] that (2.2) converges in fact for all large x 
with arg(rr) G (-7r/2, tt/2). 

3.3. Uniform representation for large x. We are now interested in the region 
of validity of (1.7) and deriving the behavior of x p (C) for large C G C. While the 
results in the previous sections used specific properties of (1.3) (positivity of hk), the 
technique of transasymptotic matching introduced in §1.2 is expected to work for 
more general nonlinear systems. We will only prove (1.7) to leading order, sufficient 
to obtain (2.7); cf. also Remark 12. The expansion (1.7) is studied in a region of x 
containing points near x^ (see (1.9)) and extending to infinity. Let e > be small, 
M G R+ and s = s(x;C) = Yle x x l l 2 jC. Define K = {x : \S(x)\ < M\s - 1| > 
e and U(s) > 1 - e}. 

Lemma 10. Let C 1 G 1R + be large. For any C with \C\ > C\ we have 
(3.14) h{x] C) = -^—+^(x-,C) 

{s-iy 

where \ip(x; C)\ < K\ lnCj" 1 in Clos(1Z), for some K not depending on C . 
Step 1: Proof of (3.14) for x G R + . From (2.2) and Lemma 7 we obtain 



(3.15) h(x; C) = h + J2 l2ks ' k + J2 12ks ~ k t 1 + r k(x)x' 1 ) k - 1 

k=l k=l 

Us + W. 12S^W s . t + £ m [(l + rt{x)x _ lf l} + 0(x _ 1} = 

^ ' k=l [x\+l 

7 7^ + v^TT^ Wlth l^i(^)l < const, (x G R fl K) 

(s — \) z in(O) 

The estimates in (3.15) are obtained in the following way: by Lemma 7 we have 
\r k (x)\ < 9/4; thus, |c fc (x)| = |(1 + r k /x) k - 1| < (e 9 / 4 - l)kx^ for k < x; the 
tail i s exponentially small, while for x G 1Z we have h = 0(x~ 4 ) and 

|x/lnC| < 2 as C — > oo. Similarly, 
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(3.16) 



±(h(x,C)-l2s(s-iy 2 ) 



const „ m , m . 

< —— for xeR + nn) 
ln(C) 



S'tep 2. Extension to 1Z. Let tp(x) = h(x; C) — 12s/ (s — l) 2 . It is convenient to change 
variables to t> = x + | lnx — ln(C/12) and define 1Z V = {v\x e 72.}. Then -0 = VO 
satisfies the system 

(3.17) m=(i, ivtvr,.?^-^ 



where = (10e" + e 2v + l)(e v - 1)~ 2 and 

(x + 2^ 2 - -e 3 " + 4e 2w + 24e"(2e" + 1) c 4 



' " 2xV2 (x + i) 2 4x 2 (e^-l) 4 x 3 / 2 x 1 /2( e ^-l)4 ^7/2 

Let [/ (£; f ) be the fundamental matrix of U' = AU (the solution with U (t; t) = I, the 
identity matrix). 2 

Remark 11. U(t;v) and £/ _1 (t;t>) are uniformly bounded on 

K M ;t = {(t,v) e TZ 2 V \ \t-v\<M} 

Proof. U is uniformly bounded since F is uniformly continuous on 1Z V . Since 
Tr(A) = we have det(U) = const = 1, so U^ 1 is bounded as well. □ 

Let v G TZ V nl. By (3.15) and (3.16) we have |^(u )| < const./ InC. We write 

(3.17) in integral form 

(3.18) iP(v) = U(v ;v)iP(v ) + U(v ;v) U~\v Q - t)R(t)dt 

The integral equation (3.18) is contractive in L°°(TIm;v ) in a ball of radius Kj InC 
for some if which is, by (3.15) (3.16) and Remark 11, independent of vq (note that 
all terms on the r.h.s. of (3.18) are in fact 0(1/ In C) if \\tjj\\ = 1). Lemma 10 is 
proven. □ 

Step 3: Pole location for large C G C. Proof of Proposition 1, (v). Let 

K = {x : \s — 1 1 = e and \v — 2Nni\ < 2n} 

If e is small enough, then for large C the curve /C is a simple closed curve in Clos(72). 
We have by Lemma 10 § K ( f h(-] C)) = — 127T2 + o(l) thus /i is not analytic in the 
interior of /C. The only singularities of h(x; C) for x ^ are double poles, cf. the 
appendix. □ 



incidentally, the equation U' = AU can be solved explicitly: for the associated second order 
equation the general solution is ip = Cit(t + l)(t - 1)~ 3 + C 2 (t~ 1 P 5 + 60tv(t + l))(i - 1)~ 3 , where 
t = e v and P5 is a quintic polynomial in t. 
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Remark 12. Using the asymptotic behavior of the solution of U' = AU and more 
careful estimates, it is possible to prove contractivity and thus (3.14) ^ n a sector of 
opening less than 7i centered on IR + (the opening cannot equals because of (2.8)). The 
asymptotic expansion (3.14) can be determined then to higher orders in a straight- 
forward way. (The equation for Gk is (3.17), with r.h.s. a bilinear expression in 
G m , m < k (cf. also footnote 2)). The result is 

(3.19) g c ~ G (s) + -G^s) + ... + ^G n (s) + ... 

with 



12s 15s 3 + 175s 2 + 30s - 2 

(3 ' 20) G ° " Gl " io a ( a -i)3 ' ' ' ' 

It is also worth noting that for all Gk(t) to be rational functions, the value —392/625 
in (1.3) is crucial. A different numerical value produces a valid expansion (3.19) but 
with Gk, k > 5 not rational. 

Remark 13. The arguments that were used to determine the asymptotic position of 
poles work even if C is not large, as long as \x\ is large, yielding the location of poles 
in the neighborhood of the imaginary axis, for large x (while on the positive real line 
C small implies x p (C) small, by Proposition 1). 

Appendix A. All solutions of PI are meromorphic. 

The fact that the only singularities of solutions of PI are double poles was indicated 
by Painleve, but it was not until the advent of the isomonodromic approach (see [1]) 
that a rigorous proof was given. For completeness, we present an independent proof 
of this crucial property of PI (and thus of (1.3), for x ^ 0), based on the idea that 
near poles, the equation PI is meromorphically equivalent to an elliptic equation [10]. 

Notations: z denotes a point where initial conditions for y are given, z will denote 
singular points of y. Let D be a disk without a discrete set of points. A(zi,z 2 ) will 
denote smooth curves in D, with the endpoints Z\ and z 2 G Clos(D) removed, such 
that the length of any subsegment A(z 3 , zA of A(zi, z 2 ) is less than 2\z 3 — z A \. 

Lemma 14. In a neighborhood of any point z e C there exists a one-parameter 
family of solutions h of (1.1), having a double pole at z, of the form h(z) = (z — 
z)~ 2 A(z) where A{z) is analytic at z, and A{z) = 1. 

Proof. Substituting y = Y^k=-2 c k( z ~ z) k in (1.1) we get 

k+2 

(A.l) (k + l)(k + 2)c k+2 = 6 ^ C mCk~m + z5 k = + 4=1 

m=-2 

so that C-2 = or 1; choosing c_2 = 1 we have c_i = and 
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k 

(A.2) {(k + l)(k + 2) - 12]c k+2 = 6 

+ z5 k=0 + 4 = i 

m=0 

Carrying on the calculation we get 

(A. 3) (c_ 2 , c_i, c , ci, c 2 , c 3 , c 4 ) = (l, 0, 0, 0, — z, c 4 ) 

i.e., c 4 is a free constant (Note: The specific form of the equation (1.1) is crucial in 
obtaining an integer power series solution. Generic perturbations of (1.1) lead to an 
inconsistent equation for c 4 . 

Finally, for some p > and all k > we have \c k \ < |(/c + l)p k+2 . Indeed, let p be 
such that |q| < |(i + l)p l+2 for i < 2. Then, assuming this property for c , Ci, c k+ i 
we have 

l (fe + l)(fe + 2)(fe + 3) fc+4 fc+4 
|Cfc+2 ' " 3 (fc-2)(fc + 5) P ~ {k + 3)P 

□ 

Lemma 15. If y is a solution of (1.1), analytic at z , then the radius of analyticity 
atz is at least R= 1/ max{|y(^ )| 1/2 , \y'(z )/2\ 1 / 3 , |yg + V6| 1/4 }- 

Proof. With y = YlT=o c k( z ~~ z o) k we have 

k 

(A.4) (A; + i)(A; + 2)c fc+2 = 6^ + z 5 k =o + S k=1 

m=0 

By the definition of R we have \cj\ < (j + 2)R~ j ~ 2 for j < 2. Inductively, if the 
inequality holds Vj < k + 1 we have 

k+2 | < <t^±m±M^. = (t + 3 )iT" 



(fc + 1)(A: + 2) 



□ 



Lemma 16. Let y be a solution of (1.1), analytic at any z e A(z ,z) and that 

y{ z ) = ^2k=-2 a k( z ~~ z ) k + 0(z — z) 5 as z — > 5 a/ong A(^o, 5), wrai/i a_2 7^ 0. 
Then y is meromorphic at z, and y is uniquely determined by z and a 4 . 

Proof. Since ^2 k= -2 a k( z ~~ z ) k * s an 0( z ~ z ) 5 formal solution of PI, for i < 3 
we must have a k = c k defined in (A. 3) (note: since y solves (1.1), the asymptotic 
expansion is differentiable). Let yo(z) be the solution given in Lemma 14 with c 4 = a 4 . 
Let v(z) = y (z) - (z - z)~ 2 and f(z) = y(z) - y (z). Then v = 0((z - 5) 2 ), / = 
0(z - z) 5 and 

( A -5) /" " (7^2 / = + 6 / 2 

or 
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f(z) = C 1 (z- ~z)- 3 + C 2 (z - ~zf - -— ^— f\s - ~z) A g{s)ds 

\ Z Z ) Jz 

(A.6) j\s - z)- s g( S )ds = Af(f) 

where g = 12fv+6f 2 = O(z-z) 7 near z. Thus C\ = C2 = 0. Let z\ G A and S be the 
Banach space of continuous functions on A(z±, z) with the norm sup A \{z — z)~ 5 f(z)\. 
It is immediate to check that if z\ — z is small enough J\f is contractive in the unit 
ball B C S. The solution of (A.6) is thus unique in B, f — 0. □ 

Let y be a solution of (1.1), and D an open disk where y is meromorphic. By 
definition, the singularities of y in D are isolated. Let z G dD be a singular point of 
y. We will show that z is a pole: £/ws implies that y is meromorphic on C. 

Let D he D without the poles of y. 

Remark 17. (i) As z — > z along A(z, z) we have \y{z)\ + | — > 00. (it) y is 
unbounded as z — > z. 

Proof, (i) Since y is singular at z then by Lemma 15 y(z) or |/(z) must be large 
for z near z. (ii) Multiplication of (1.1) by y' and integration yields: 

(A.7) {y'f = % 3 + 2zy - 2 f y(s)rfs + C 

J a 

If y were bounded along A, then by (A.7) y' would be bounded as well. □ . 

Proposition 18. Let z G D be a regular point of y and A = A(z , z) C D. We have 
y~ 3 (s) f^y(t)dt — > 0, \y(s) \ — > 00, and \y'(s)\ — > 00 as s — > 5 on A. 

Proof. For any continuous function / unbounded at 5 we have 

liminf |/|~ 3 / \f\(t)d\t\ < liminf |/|" 3 max |/| length(A) = 

z^z J a z->z A(a,z) 

Thus, if we had y~ 3 (s) f*y(t)dt 0, then for any 7 > small enough 
(A.8) I ^ y{t)dt\= 1 \y\s n )\ 

J a 

for some sequence s n — > z on A. If 7 < |, by (A.8), (A.7) and Remark 17 we have 

(A. 9) \y(s n )\—>oo and \y\s n )\ -> 00 and / — > 00 

Consider now the solution ft, of the equation 



(A.10) 



(ft') 2 = 4h 3 + 2s n h + C 
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satisfying h(0) = y(s n ),h'(0) = y'(s n ) (with C chosen accordingly in (A. 10)). It 
follows from (A. 9) that C = C — 2 J^ n y(s)ds — > oo. Define the function H by 

aH{(3z) := h(z), where a = C 1/3 , (3 = C 1/6 

It follows that 

(A. 11) (H'f = AH 3 - 26H + 1, where 9 = C~ 2/3 

H is thus an elliptic function with with periods n 12 = II + o{6) for some II E C 
(independent of 0). By (A.8) \y(s n )\ = (2-y _1/3 ) | C' 1/3 1 ( 1 + o(l)), then using (A.7), we 
get \y'(s n )\ = (4 7 ~ 1 e^ + 2)- 1 / 2 1 C' 1 / 2 1 (1 + o(l)) for some 0„ E R. 

Therefore \H(0)\ = (27)- 1 / 3 + o(l), and 1^(0)1 = (47-V^ + 2)" 1 / 2 ! + (1). With 
the substitution aY(/3z) := y(s n + z) we have 

{Y'f = AY 3 + 29Y + 1 + 2e - jf Y(t)rf^ , where e = C^ 5/6 , 
(A.12) or Y" = QY 2 + 6 + es 

We substitute Y = H o G; G(z) = z + g(z) in (A.12) with g(0) = g'(0) = and 

2g'H" oG+ (g') 2 H" o G + g"H' o G = ez 
(A. 13) or 2g'(H' o G)' + g"H' o G = ez + (g') 2 H" o G 

which can be rewritten in the integral form (g(z) = f A ^ Q _^ w(t)dt, G = z + g): 
w = (H'o G)- 2 [ H'{G{u)) [eu + H" \G(u))w 2 (u) 2 ] du = Af(w) 

JA(0,z) 

(A.14) 

Let A e be the disk centered at 0, of diameter jej" 1 / 4 , and A e 3 be A e without 
the disks of diameter ^(7) around the poles and zeros of H'{z). Let A C A e , let 
a E (1/4, 1/2). Let S be the Banach space of continuous functions on A in the sup 
norm and B = {w : ||w|| < |e| a } C S. If w G B, then \g\ < \ f A{0 z) w(t)dt\ < \e\ a ' 1 / 4 . 
B is preserved by N provided that 

(A.15) |e| Const(-f)\e\- 1/2 + Const( 7 )|e| 2a " 1/4 < |e| a 

which holds for small e. Similarly for u,v E B, \\Af(u) — Af(v)\\ = 0(e a_1 / 4 ) \\u — v\\. 
N has thus a unique fixed point in B. 

Therefore \g(z)\ < |e| a_1 / 4 for z E A e and thus (when e is small) 

(A.16) sup \a- l y{s n + /T 1 -) - H{-)\ = sup \H{G{-)) - H{-)\ = 0(e a ~ l l A ) 
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Note now that as s n — > z we have f3~ l = C^ 1 ^ — > while a — > oo and | / 9 _1 e~ 1 / 4 | oc 
C* 1/24 -> oo. Thus L = s n + a _1 A e is an a: -1 -lattice for D (i.e., Vz <E D , 3zi G L : 
|z — zi\ < and (A. 16) implies y is unbounded near zq, contradiction. 

Since y~ 3 (s) f s y(t)dt — > 0, the fact that min{]y|, \y'\} — > oo now follows from 
Remark 17 and (A. 7). □ 

Proposition 19. z is a second order pole. 

Proof. The substitution y = 1/v 2 transforms (A. 7) into 

(A.17) {v'f = 1 + \v{zfz - \v{zf f v(sY 2 ds + Cv(zf 

2 2 J zo 

It follows from Proposition 18 that along a curve A(z , z) avoiding the poles of y we 
have v 6 J* o v- 2 (s)ds -> 0. With ( = z - 5, (A.17) implies (w') 2 = 1 + o(l) so that 

= ( 2 (l + o(l)) 2 near 5. Using (A.17) again we now get v 2 = (C + EL2 ^fcC fe + o(C 6 ))- 
Lemma 15 implies that y is meromorphic at z. □ . 

Acknowledgments. The author would like to thank Professors Percy Deift and 
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